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lf~) ■ We consider radiative corrections to the decay rate F(H — > ZZ) of the heavy 
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q^ ■ perform a one loop Feynman diagram calculation in the on-mass-shell renormalization 



C-P-even Higgs boson of the minimal supersymmetric model to two Z bosons. We 



scheme, and include the third generation of quarks and squarks. The tree level rate 



j^. is suppressed by a mixing angle factor and decreases as 1/Mjj for large Mh- The 



corrected rate overcomes this suppression and increases with Mjj for Mj/>500 GeV. 



S^ The corrections can be very large and depend in detail on the top squark masses and 

^4-term, as well as the supersymmetric Higgs mass parameter \i. 
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I. INTRODUCTION 

One of the least attractive features of the standard model (SM) is the existence of the 
naturalness problem. Roughly speaking this means that when one computes corrections to 
the Higgs boson mass one finds quadratically divergent contributions. This situation implies 
that input parameters must be extremely fine-tuned at high energies to yield the low energy 
physics that we observe, a situation that is unappealing especially in connection with GUTs. 

One way to control the naturalness problem is to consider supersymmetric (SUSY) exten- 
sions of the standard model. Here the quadratic divergences are cancelled by loop diagrams 
involving the superpartners of the SM particles. We know that SUSY must be broken in 
the real world, and yet the scale of supersymmetry breaking must not be too large or the 
hierarchy problem will be reintroduced. Thus, although superparticles must be sufficiently 
heavy to have avoided detection at present colliders, they cannot be much heavier than a 
few TeV if we are to meet the naturalness criterion. 

In this work we will be concerned with the simplest supersymmetric extension of the 
standard model (MSSM) |],[|. In the MSSM we need two Higgs doublets H ± and H 2 to give 
masses to up and down type fermions and to assure cancellation of anomalies. The neutral 
Higgs spectrum consists of two CP-even Higgs scalar particles H and h (where Mh > M^), 
one CP-odd particle A, and a Goldstone boson G which is "eaten by" and gives mass to 
the Z boson. The Higgs sector of the MSSM is highly constrained. At tree level the Higgs 
boson masses and couplings are determined by two input parameters. We take these to be 
the mass of the CP-odd Higgs boson Ma and an angle (3 which at tree level is given by 
tan/? = 1)2/ V\ where v 2 and V\ are the vacuum expectation values of the two Higgs boson 
fields H 2 and H\. The tree level masses of the CP-even Higgs bosons are then given by 



M 2 H , h = 1 -[m\ + M 2 z ± ^{Ml + MIY - 4M1M1 co S 2(2/3)) . (1) 

The above equation implies the inequalities Mh < Mz, Mh > Mz and the sum rule 

M H + M% = M z + M\. 



Recently it was shown that one loop corrections involving top-quark and squark loops 
can significantly modify the sum rule || and also violate the bound M^ < Mz @|J- For 
1 TeV squark masses the correction to the light Higgs boson mass is of the order 20 (50) GeV 
for a top mass of 150 (200) GeV. Corrections to the neutral Higgs boson mass sum rule due 
to the gauge-Higgs and gaugino-higgsino sectors were considered earlier and were found 
to be generically small. 

In this work we consider corrections to the decay rate T(H — > ZZ) which is relevant 
for the detection of the heavy Higgs boson at a proton supercollider such as the SSC via 
the "gold-plated" mode H — > ZZ — > £ + £~£ + £~, where £ is e or fi. We confine ourselves 
to corrections due to third family (top and bottom) quark and squark loops. Previous 
work on this subject has appeared in Ref.'s |7[ and || where the effective potential and 
the renormalization group methods are used. We perform a Feynman diagram calculation 
utilizing the on-mass-shell renormalization scheme, and present explicit analytic results. The 
structure of the paper is as follows: in Section 2 we present our renormalization procedure, in 
Section 3 we discuss our results, Section 4 lists briefly our conclusions, and in the Appendix 
we present the necessary explicit formulas. 

II. FORMALISM FOR RADIATIVE CORRECTIONS 

Due to the presence of mixing in the CP-even and CP-odd sectors the renormalization of 
the Higgs sector of the MSSM presents a few complications when compared to the standard 
model. Therefore, in this section we present in detail our renormalization procedure. We 
follow the approach of Aoki et. al. [[J adapted to the MSSM. 

The Higgs potential in the MSSM is 



V = - +9 (h[*H[ - H l 2 *H l 2 ) 2 + —{H^Hil 2 (2) 

g \ '2 



+ {m\ + n 2 )Hl*H[ + {m 2 2 + [i 2 )E^E\ - (m^tfji^' + h.c. 

where g(g') is the SU(2)l(U(1)y) gauge coupling, the m;'s, (i = 1,2,3) are the soft super- 
symmetry breaking Higgs sector mass parameters, and /i is the supersymmetric Higgs mass 



parameter. We can absorb fi 2 in Eq.(l) by redefining m 2 + fi 2 — > m^ and similarly for m 2 . 
Hi and H 2 are given in terms of the shifted (but unrotated) fields by 

1 (v 1 + S 1 -iP 1 \ 1 ( V2H, 

V2\ ^2H_ )' 2 V2\v 2 + S 2 + iP 2 , 

In order to discuss the tadpole and mixing structure of the theory we need the terms that 
are linear and quadratic in S±, S 2 and quadratic in Pi, P 2 . These are given by 

\r I 9 "r 9 ( 2 2\ , 2 2 \ a i I 9 ' 9 I 2 2\ , 2 2 \ o 

V s = I [Vi - v 2 )vi + m x vi - m 3 v 2 1 Si + I (u 2 - ^ ju 2 + m 2 v 2 - m 3 t;i I 6 2 

(g 2 + g' 2 ( ~ 2 2n . ™?\ c2 . (g 2 + g' 2 f o 2 2 n . ^ c2 

+ (^-(^i " « 2 ) + -yj *i + (^-(^ - vi) + —J S 2 

2 j. „/2 \ 

^1^2 + ml SiS 2 (3a) 



0+0' .. .. , .. 2 



V 4 
v„ = (^r^K - -I) + "y ) ^i + (^r^ - *?) + t) Pl ~ mlPlP2 - (3b) 

We now define the coefficients of Si and S 2 in Eq.(3a) to be 

Ti = ^-^(v 2 - vl)vi + mjvi - m\v 2 (4a) 

T 2 = ^(v 2 - v 2 )v 2 + m^ 2 - m\vi. (4b) 

Eliminating to^, m 2 in favor of Ti, T 2 from Eqs.(4) and substituting back in Eqs.(3) we 
obtain, using a matrix notation 

V.-(ft*)g) + i(ft*)(; J) (*) (-) 

1 / ^«? + ^ "^-^\ /* 

+ 2 ^ ^ \ -^v lV2 - m^ £t£„2 + ro§£ A ^2 . 

^ \Po 2 V -m* m^ / V p 2 



The next step is to introduce rotation matrices 0(a) and 0(/3) such that the part of the CP- 

even and CP-odd mass matrices that does not depend on Ti, T 2 is diagonalized. Specifically, 

by defining 

Si \ ( H \ ( cos a — sin a \ ( H 

= 0(a) I } = [ 
S 2 J \ h ) \ sin a cos a J \ h 



and 



we find that 



PA fG\ /cos/? -sin/A /<T 

P 2 J \AJ V sm /3 cos P J \ A 



v,= (Hh) ^y- {Hh) o ( - a) ^ ^ 0{a) Q + \ (Hh) ( M « ;j(^ 



(6a) 



V p ^(GA )0( - m ■ t JO W[aJ+ - ( GA )[o U |. (Ob) 



1 (^ o \_ _ /g\ | /o o \ /<r 

2' \ n Ik / \ j / 2 V \ n A/T2 

•(■7 



Here we have defined 

TA (T H 

= 0(a) 

The parameters /3, a, M#, M^ and M A are related to the original fundamental parameters 
Vi, V2 and m§ by the following formulas 

tan/?=— , M\ = ml {tan 13 + cot 13), tan 2a = ^ + -f tan 2/3, (7) 

"Uj -IV-l a — IVl y 

as well as Eq.(l). Here we used M\ = 9 + ® (v\ + f|). Carrying out the remaining matrix 
multiplications involving the tadpole contributions to the mass matrices we obtain the final 
result 

1 (M 2 H + b HH b Hh \ (H\ 

V s = HT H + hT h + -(Hh) (8a) 

2 V b Hh M 2 h +b hh ) \h) 

1 (b GG b GA \ (G\ 

V P = -(GA)[ (8b) 

2 \b GA M\ + b AA ) \A) 



with 



bnH = t. [Th(cos 3 a sin (3 + sin 3 a cos (3) +T h sin a cos a sin(o; — (3)) 

v sin2p v 

sin 2o; 
fefffc = — ^tt. (T H sin(a - /3) + T h cos(a - (3)) (8c) 

2 

&Wi = 77 (Th cos a sin a cos(o; — f3) +T ft (cos 3 a cos (3 — sin 3 a sin (3) 

v sin 2p 



and 

bcc = ~ (T H cos(a - 0) - T h sin(a - /?)) 

v 

bo a = - {T H sin(a - (3)+T h cos(a - (3)) (8d) 

b>AA = 7. (Tff(sin 3 /3coso! + cos 3 /3 sin a) + T/ l (cos 3 / 9coso; — sin 3 /? sin a 

v sin2p v 

The terms linear in H and h are to be thought of as counterterms for the tadpoles. To each 
order in the loop expansion we require that the total tadpole contribution vanishes. At tree 
level this implies —iTn = = —iTh- This then gives the conventional tree level masses. At 
one loop —iT H (—iT h ) must cancel the one loop H (h) tadpole diagrams ir H (ir h ) (Fig.l). 
These conditions determine Th and T^ and Eqs.(8) determine their contribution to the one 
loop mass matrices. 

Taking as renormalized inputs tan/3 and Ma we calculate the physical masses Mh, M^ 
and the decay rate T(H — > ZZ) at one loop. It follows that the measurement of any two of 
the physical quantities Ma, Mh, M^ and T(H — ► ZZ) will allow us to make a prediction for 
the other two. We stress that (3 is only to be viewed as a useful parametrization of physical 
observables. Since by itself f3 has no physical meaning we can renormalize it in any suitably 
convenient way. We explain our renormalization prescription for /3 below. 

From this point on we adopt the following notation conventions: a quantity such as a 
field, coupling, or mass with a subscript indicates a bare quantity, renormalized quantities 
have a subscript r, and physical observables such as the pole of a propagator do not have 
subscripts. The bare tree Lagrangian contains 

C D ^Hod^Ho + ^hod^ho (9) 

-\{M 2 Ho + b HH )Hl - l -{M 2 ho + b hh )hl - b Hh H h 

where M Ho and M| o are taken to be functions of M^ , (3 and M Zo as given by equation 
(1). We now write the bare parameters in terms of renormalized parameters and shifts 

[3 = f3 r + 5[3, M 2 Ao = M 2 Ar +SM 2 A , M 2 Zo = M 2 Zr + 5M 2 Z (10) 



and also introduce wave function renormalization 
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Ho = Zjj H H r + Z?r h h r} h = Z? h h r + Z? H H r . 



(11) 



Note that Zj 1H = 1 + 0(a), Z£ h = 1 + 0(a) while Z Bh , Z£ H , b H H, b H h, and b hh are all 
0(a). Substituting equations (10) and (11) into (9) we obtain the one loop renormalized 
two-point functions 

i dM 2 

iThh(p 2 ) = (Z HH ) 2 (p 2 - M 2 H ) - ^^Sxi - b HH + Tl HH (p 2 ) 



i dM 2 

iT hh (p 2 ) = {Zl h ) 2 {p 2 - M 2 hr ) - -^Sxi - b hh + U hh (p 2 ) 

1 1 

,2\ 72 I Ji S/f2 \ i 72 I Ji S/f2 \ L i TT { J2\ 



(12) 



zT Hh (p 2 ) = Z Hh (p 2 - M Hr ) + Z'UP 2 ~ Ml) - b Hh + U Hh (p 2 

where X{ r = M\ , M§ r , /3 r and the IPs are the scalar self-energies (Fig. 2). The on-shell 
renormalization conditions are [01 



iT HH {M H ) = iT Hh {M 2 ) = iT Hh {M 2 H ) = iT hh {M 2 ) = 



dY 



HH 



dp 2 



l = i 



or 



iii, 



p 2 =Mf 



dp 2 



(13) 



p 2=M l 



Here Mh and Mh are the physical masses of H and h. Making the definitions 5M H 
Uhh(M h ) — buH and similarly for 5M%, we obtain from Eqs.(12) and (13) 



dM 
M 2 = Ml + „ 
H Hr dXi 



H r 



5xi - SM H 



(14a) 



, > , r> 9M h 



Ml = Ml + ~-rr^5x t - 8M 2 



(14b) 



z HH = i- \n' HH {M Hr 



4 = i '"' " r2 



;K h (M 2 h 
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Ml -M 2 



-n Hh (M 2 ) + b Hh 



(14c) 
(14d) 
(14b) 



Z HH = M2 _ M2 (-n^(M^) + b Hh ) , (14f) 

fir Hr 

where the prime in Eqs.(14c,d) indicates differentiation with respect to p 2 . Note that M]j 
and M| have the same functional form as in Eq.(l) except that they are functions of 
renormalized quantities, i.e. 



M 2 Hr , hr = g [ M l + M l ± V( M i + M i) 2 - 4M l M i cos 2 (2/3 r ) J . (15) 

We now drop the subscript r on M^ r , M^ and (3 r . Eqs.(14a,b) determine the physical 
CP-even Higgs boson masses in terms of self energies, tadpole contributions, and shifts of 
the inputs parameters 5x{. We now determine the shifts. The shift 5M A is defined so that 
M A is equal to the physical A mass. An analysis similar to that of the CP-even sector yields 

5M\ = U AA (M 2 A ) - b AA . (16) 

Additionally, we find for the shift in the Z-boson mass 

5M 2 Z = lf zz {Ml) (17) 

where n^ z is the transverse part of the Z boson self energy, 11^ = g^U^z + ^-H zz . At 
this point it is worth noting that if we are only interested in the sum M\ + M\ we do not 
need a specification for 5(3 . When Eqs.(14a) and (14b) are added the terms proportional 
to 5(3 cancel leaving 

M 2 H + M 2 h = M\ + M 2 Z -U HH (M 2 H ) - U hh (M 2 h ) + U AA (M 2 A ) + U T ZZ (M 2 Z ) 

+b HH + b hh - b AA (18) 

This is just the renormalization of the neutral Higgs boson mass sum rule and the divergences 
in Eq.(18) implicit in the Il's and 6's cancel leaving behind a finite correction. Since we 
demand that M# and M^ are physical masses they must be individually finite. Equivalently, 
since Mfj + M\ is finite we must have that M^ — M 2 is also free of divergences. This latter 
requirement gives 

7 



<9A <9A <9A 

W 5p + mt Ml + m\ ml ~ ml + ml = finite (19) 



where A = J(M\ + M§) 2 — 4M\M Z cos 2 (2(3) . The above equation clearly determines only 
the "infinite" part of 5(3. By "infinite" we mean the part that is proportional to Cuv — 
- — 7 + log Att in dimensional regularization. To fully specify 5(3 we take a MS- type approach 
and define 5(3 to be purely "infinite" so that Eq.(19) becomes 

W 1 " - ~ (m ml + ik m ^ ~ SM » + SM! ) . (20) 

where the subscript oo on a quantity indicates the "infinite" part of that quantity. Eq.(20) 
implies 



5(3 = = — = — - (21) 



1 

2MiMfsin(4/3) 

(M\ 5 Ml + Ml 5 M 2 Z ) cos 2 (2(3) - M 2 H 5M 2 h - Ml 5M 2 H 



x 



This definition of (3 at one loop gives renormalized CP-even Higgs boson masses in close 
agreement with those obtained using the effective potential J3J. This shift in 5(3 induces a 
shift in a through equation (7) 

■ ,* J &P M 2 A 5M 2 z -M 2 z 5Ml\ 

We now come to the renormalization of the HZZ coupling. The bare HZZ and hZZ 
couplings are given by 

^ ZZ = A# '° 2 M| ° 2 ri cos(/? - ao), \f z = \» zz tan(/? - a ). 
M Wo (M 2 Zo - M^ o )2 

Defining 

i ii 

e = e r + 5e, Zfj H = 1 + 8Z\ H , {Z zz ) 2 = 1 + <5Z ZZ 

i i 

(here Z^ = Z'^Z^ + Z ZA A^ r where Z^ (Z^ r ) is the bare (renormalized) Z boson field 

and A^ r is the renormalized photon field) we obtain for the renormalized one loop 3-point 

function 

? H J Z = (*r™ + AgP)fc* + ATjf* (23) 



where \? zz = ^cc^-o,) and 

M Wr (M 2 Zr -M^ r )^ 

,hzz .HzzfSe^SMl 1 8MZ-8M& 1 6M& 

AcT ~ K U + 2 Mf 2 M z -M^ 2 Ml [ ] 

i i 

— tan(/3 r — a r ) (5/3 — 5a) + 5Zfj H + 5Z ZZ + -%? tan(/3 r — a r 

and Ar^ zz is the explicit one loop Feynman diagram contribution (Fig. 3). The angle 

a r is defined as in Eq.(7), but with the right hand side written in terms of renormalized 

i i 

quantities. The expressions for 5M Z , 5/3, 5a, Zfj H and Z£ H in terms of self energies and 

tadpole contributions are given in Eqs.(17), (21), (22) and (14c, e). We simply state the 

results for the remaining shifts 5e, 5M^ and 5Z ZZ . We have 

^4<(o )+ (^m (2B) 



e 2 77 v ' \ As w c w J M. 
5M 2 W = n^(M^), 5Z ZZ = -Yi T zz \M\ 



where c w = M w /M z and s w = J 1 — c%,. We note that 11^,(0) vanishes in our case. The H-h 

i 
mixing gives a contribution to T^ zz through the term proportional to Z£ H . The quantity on 

the R.H.S. of Eq.(23) is given as a sum of terms which are individually divergent. In the full 
sum the divergences must of course cancel. We checked both analytically and numerically 
that this is indeed the case. The renormalizability of the theory requires that the definition 
of 5f3 which renders the CP-even Higgs boson masses finite also gives finite couplings. 

The explicit one loop Feynman diagrams shown in Fig. 3 give a contribution to the three- 
point function which can be expanded in terms of form factors as 

AT% ZZ = D g^ + D lP ^ 1 + D 2 p%p» + D 3 rfp" 2 + D^p\ (26) 

(a form factor proportional to e tlua,3 piaP2f3 vanishes by CP invariance). The formula for the 
decay rate at one loop is 



r= 12 ^^ {(l-4r + 12r 2 )((Af^) 2 + 2A r ™Re(AgP + D ) + \\% z + D f 
+ M 2 H (1 - 2r)(l - Ar) (\* zz Re(L> 4 ) + Re[(\^ zz + D )D* 4 ] 
+ MJ { ^--2t) 2 \dA (27) 



where r = M^/Mfj and we list Do and -D4 in the Appendix. 

We note that the terms in the above expression which do not involve X^ zz are formally 
of 0(g 6 ). Nevertheless we find that for large Higgs boson mass (Mg 3> Mz) they are numer- 
ically important. This is because \f? zz is proportional to cos(a — (3) which is proportional 
to \jM\ for large M# and hence small. Keeping these 0(g 6 ) terms is consistent: the terms 
in the amplitude that are of 0(g 5 ) which arise at two loop level also give a contribution of 
0(g 6 ) in the decay rate, but these two loop 0(g 6 ) terms are proportional to cos(a — (3) and 
are thus suppressed when Mh 3> Mz, in precisely the region where the 0(g 6 ) terms in our 
one loop expression become large. 

III. RESULTS 

In the MSSM at tree level the decay rate T(H — > ZZ) is suppressed relative to the same 
decay rate in the standard model by the factor cos 2 (a — 0). The "gold-plated" decay mode 
H — > ZZ — ► 4£ has great discovery potential for a standard model Higgs boson at a proton 
super collider such as the SSC for Higgs boson masses 130 GeV < M^ < 800 GeV |ID| . The 
discovery potential for the heavy Higgs boson of the MSSM in this mode is not as promising 
due to the above mentioned suppression factor. However, the "gold-plated" mode may be 
the only discovery mode for the heavy Higgs boson at a hadron collider ||11| . The discovery 
potential is improved when radiative corrections are taken into account. 

We discuss our numerical results below. We have checked our numerics in a number of 
ways. First, we checked the cancellation of divergences as mentioned in the last section. 
Second, we found our result for the correction to the neutral Higgs boson mass sum rule 
agreed very closely with that of Ref. 0|. Lastly, we checked that our calculation, when 
modified to give the correction to the standard model Higgs boson decay rate to two Z's 
due to an extra heavy fermion doublet, agrees with the results of Ref. [[14] . 



In Fig. 4a we show the tree level and radiatively corrected decay rate versus the heavy 
Higgs mass for tan /3=5 and a top quark mass of 160 GeV. In this figure we have not included 
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mixing effects, i.e. A t — Aj, — \x — and the squark masses are all equal. We show the 
corrected rate for the two squark mass choices M sq = 300 GeV and M sq = 1000 GeV. We see 
in Fig.4a the importance of keeping corrections which are of 0(g 6 ) in the rate. The one loop 
corrections which contribute 0(g 4 ) to the rate fall with M# (as they multiply the tree level 
coupling). However, the one loop corrections which contribute 0(g 6 ) to the rate increase 
as Mu increases. Hence, these terms eventually dominate the rate as Mu becomes large. 
In Fig. 4a the corrected rate is dominated by the 0(g 4 ) terms for small M#, and hence it 
initially falls as Mu increases beyond the kinematic suppression. Eventually, however, the 
terms of order 0(g % ) become larger than the 0(g 4 ) terms and the rate then rises with Mu- 
This begins to occur for values of Mu of about 500 GeV. 

In Fig. 4b we show the rate versus tan (3 for a Higgs boson mass of 300 GeV, a top quark 
mass of 160 GeV, a squark mass of 1 TeV, and again for no mixing. We see that the corrected 
rate is approximately twice as large as the tree level value, almost independent of tan/3. As 
we will discuss below, the rate depends dramatically on tan/5 once mixing is included. 

In Fig. 5 the ratio of the radiatively corrected rate to the tree level rate is shown versus 
the top quark mass, for the same set of parameters as Fig. 4b, and tan/?=5. Fig. 5 illustrates 
that the corrected rate depends strongly on two parameters in the case of no mixing. Clearly 
the rate depends on the value of the top quark mass. But note for M#=l TeV that even for 
a top quark mass as small as 100 GeV the corrected rate is still over a factor of two larger 
than at tree level. Thus the relative size of the correction depends greatly on the value of 
Mu as well. Note, however, that when the top quark mass is less than around 120 GeV we 
expect that the corrections from other sectors will be of the same order of magnitude as the 
correction due to the quark/squark sector included here. 

When mixing is included the parameter space increases. We will choose a point in mixing 
space and examine the effect of mixing in deviations from that point. We choose A— terms 
A t = At, = 600 GeV and squark masses m tl = 771^=600 GeV, and fh t2 = m;, 2 =300 GeV. 
Additionally, we will consider the two cases \x = ±400 GeV. In all three of the figures 6, 
7 and 8 the heavy Higgs boson mass is set to 300 GeV and the top quark mass is 160 
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GeV. In order to isolate the effect of mixing we will plot the ratio of the corrected rate 
including mixing to the corrected rate with no mixing (where the common squark mass is 
set to 600 GeV). In Figs. 6 we plot this ratio vs. the squark mass mt x . We find that the effect 
due to mixing is strongly dependent on tan f3 and \x. For large values of tan f3 the effects of 
mixing are greatly enhanced. As shown in Figs. 6, the inclusion of mixing can change the 
rate by a factor 1.3 for tan/3=2 and for tan/3=20 by a factor 2.7 or 0.3, for /i=-400 GeV or 
/i=+400 GeV, respectively. 

Similar ratios are seen in Figs. 7, where the ratio of the corrected rate including mixing to 
the corrected rate with no mixing is shown vs. A t , the top squark mixing parameter. As in 
Figs. 6 the two curves for // = ±400 GeV are similar when tan/3=2; the rate can be increased 
by 50% or decreased by 25%. If tan/3=20 the effects of mixing are more pronounced and the 
ratio varies between roughly 1/3 and 3. The /z=400 GeV curve in Fig. 7a (and the tan/3=2 
curve in Fig. 8) does not span the entire ordinate axis shown because an unphysical region 
of the squark mixing parameter space is encountered. In Fig. 8 we plot the (mixing) to 
(no mixing) ratio vs. the supersymmetric Higgs mass parameter fx. We see there is little 
dependence on \x for small tan/3, while for larger values of tan/3 the dependence is quite 
significant. If tan/3=20 the ratio varies between 4 and 1/36 as \x varies from -750 to 750 
GeV. Finally, we note that there is very little dependence on the bottom squark masses and 
A— term A^ for the mixing configurations considered. 

IV. CONCLUSIONS 

To summarize, we have computed the one loop corrections to the decay rate T(H — »■ ZZ) 
in the MSSM including third family quark and squark loops. We perform a Feynman diagram 
calculation in the on-mass-shell renormalization scheme. As the tree level rate falls like 1/Mh 
for large Mh and we find corrections that grow with Mh, the corrected rate may be many 
times the tree level rate. For example, at Mh = 1 TeV the corrected rate may be 13 times 
the uncorrected rate for m t =200 GeV (with no squark mixing). The corrected rate depends 
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very strongly on the squark mixing parameters. For example, for the mixing configuration 
considered here, the rate varies by two orders of magnitude as the Higgs mass parameter \x 
varies between ±750 GeV. Indeed, the squark mixing parameters fi, A t , and the top squark 
masses, in addition to the top quark mass, must be measured in order to test the Higgs 
sector of the MSSM. 
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In this Appendix we give explicit analytic expressions for the self energies, tadpoles, and 
form factors introduced in the text. Our expressions are given in terms of the standard 
A, B, C functions introduced by Passarino and Veltman |T3] which appear in one loop cal- 



culations. We adopt the metric (1,-1,-1,-1), which is different than that of Ref. jT2|. Explicit 
analytic formulas for these functions appear in Ref. [pl| . 



To make the equations more concise we adopt the following conventions. N c denotes 
the number of quark colors. The index a runs over the top and bottom sectors while the 
indices i, j, and k run over squark mass eigenstates. Thus, m a denotes a quark mass 
while fh a i denotes a squark mass. For the A and B functions we define A a = A{m 2 a ) 1 
A ai = A{m 2 ai ) } B 0a = B (p 2 ,m 2 a ,m 2 a ), B 0aij = B (p 2 , m 2 ai , fh 2 aj ) and similarly for the rest 
of the .B's. AC function has six arguments: three external squared momenta and the 
three squared masses of the particles which appear in loop of the 3-point diagram. We thus 
define C 0aijh = C (Mf, Mf, M 2 H , m 2 ai , fh 2 aj ,m 2 ak ) and C 0a = C (MlMlM 2 H ,m 2 aJ m 2 a ,m 2 a ) 
with analogous definitions for the rest of the C's. 

First we give expressions for the Higgs boson self energies. 
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IW P 2 ) = n c Y,(v a %) 2 B 0aiJ -N C Y, u^a m (Al) 

aij ai 

-12N c ^(V Q H ) 2 (mlB 0a + p 2 (B 21a - B la ) + -i(m 2 - ''' 



48tt 2V -" 6 

The various V and £/ vertex factors are shown in Fig. 9 and explicit expressions appear in 
Refs. [|],|nj. However, the H — h — %l — qkL and H — h — q~kR — q~kR vertices given 



in Ref. [10] are incorrect. In the notation of Ref. [10 1 the above couplings are 



%g l sin 2a / T 3k - e k sm 2 W m\ \ ig 2 sm2a 2 m\ 

a 2 Th Tit Dk and a 2ek tan °™ ~ -rjt D 

4 V cos 2 6» w Mfa ) 4 V A% 



respectively (D up = 1/sin /3, D down = -l/cos 2 /3). 

U hh is given as U HH with V^j -> V^, U?f -> U^, and V£ -> V£. H-Hh is given as 
U HH with (V^) 2 - i^V£, £/*f - E®, and (Vf ) 2 - V^ V a \ ILu is given as U HH 
with V^- -> V^-, U£? -> C/Af, V; H -> K/, and 5 0a -* iS 0o . Next we list the transverse 
part of the gauge boson self energies. 

nL(p 2 ) = N c ^u^A m 

ai 

- 2N c Y / (V^) 2 (rhl l B 0aij - (rh 2 m - m% + p 2 )B lmj + p 2 B 21mj + 

aij 

~ 8iV c £(V/J 2 m^ + ((V a z ) 2 + {Vi a ) 2 ) p 2 (B 21a - B lo 

a ^ 

^wwiP 2 ) = N c J2uZ W A ai (A3) 

ai 

-i ~ 9 ~ 9 9 

- ^/V c 2-A v ij ) m ti B - [m ti -m bi +p )B 1 + p B 21 + — ^ 





(A2) 


1 m 2 ai + m 2 aj 


-£) 


16tt 2 ^ 2 



16^ 2 6 

^2/_2 D W t^2 _2 | o^nW i o„2 D 
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- 8JV C (0 m 2 ^ - (m 2 - m 2 + 2p 2 )B 1 M ' + 2p 2 B 



where B w = B(p 2 ,m 2 ,ml) and B w = B(p 2 ,m 2 ,ml). l£E is given as Ii T zz with V^„- 



V^jj, C/^f — > C/J2; Kf ~ * V^7, Vg| — ► 0. The heavy Higgs boson tadpole contribution is 
given by 

T H = N c Uj2 V a H m a A a - Y, V^A^ . (A4) 

^ a ai 

Th is given as T# with V^| — > V^ u and V^ — > V£- Lastly, the two three-point Feynman 
diagram form factors which are relevant for calculating the Higgs boson decay rate are 
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D = 8N c ^m a V a H ((V z ) 2 + (V z a 



Z\2\ 



M 2 
x ( 4C 24a + (M 2 H - 2M 2 z )C 12a + 2M 2 z C lla - -^C 0a - B (M 2 Z , m 2 a , m, 



SN c J2m a V a H (V 5 z a ) 2 ((M 2 H - Ami) C 0a ~ 2B (M 2 z ,m 2 a ,m, 



8iV c £ V^V* l3 V z jk C Uai]k + iV c £ V^U zz B (M 2 h , rh 2 ai , rh%) 

aijk 



a i j 



(A5) 



and 



D 4 = 8N c Y,m a V^(((V a z ) 2 + (V z a ) 2 ) (4C 23a + C 0a - 4C 12 J - 2(V z a ) 2 (C lla - C 12a ] 



- 8N C 22 V aki V aij V ajk [CWciijk ~ ^^ a 
aijk 



ijk 



(A6) 
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FIGURES 
FIG. 1. The Higgs boson tadpole diagram. 

FIG. 2. The Higgs boson and gauge boson self energy diagrams. 

FIG. 3. The one loop H — ► ZZ Feynman diagrams included in our calculation. The dashed 
loops represent squarks, the solid loop represents quarks. 

FIG. 4. The tree level and one loop corrected decay rate Y(H — > ZZ). Fig. 4a shows the rates 
vs. the heavy Higgs boson mass for tan/3 =5, m£=160 GeV, and no mixing. Fig. 4b shows the 
decay rates vs. tan/3 for the same parameters as in Fig. 4a, with Mh=300 GeV. 

FIG. 5. The ratio of the corrected decay rate to the tree level rate vs. the top quark mass. The 
parameters are the same as in Fig. 4a, with the Higgs boson mass set to 300 GeV and 1 TeV. 

FIG. 6. The ratio of the corrected rate including mixing to the corrected rate without mixing 
vs. the top squark mass m^ as explained in the text. The top quark mass is 160 GeV. 

FIG. 7. The ratio of the corrected rate including mixing to the corrected rate without mixing 
vs. the top mixing parameter At, as explained in the text. The top quark mass is 160 GeV. In 
Fig. 7a the point where the /x=400 GeV curve stops corresponds to an unphysical point in squark 
mixing parameter space. 

FIG. 8. The ratio of the corrected rate including mixing to the corrected rate without mixing, 
as explained in the text, vs. \x. The heavy Higgs boson mass is 300 GeV, and the top quark mass 
is 160 GeV. The curve for tan/? = 2 stops at an unphysical point in the squark mixing parameter 
space. 

FIG. 9. The vertices used in the Appendix are displayed. The values of the vertex factors 
may be found in Refs. |i"l,|i"of. Note that the value of U£-j listed in Ref. [10] is incorrect (see the 
Appendix for the correct values of these vertices). 



